Loose packings of frictional spheres 
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We have produced loose packings of cohesionless, frictional spheres by sequential deposition of 
highly-spherical, monodisperse particles through a fluid. By varying the properties of the fluid and 
the particles, we have identified the Stokes number {St) - rather than the buoyancy of the particles 
in the fluid - as the parameter controlling the approach to the loose packing limit. The loose packing 
limit is attained at a threshold value of St at which the kinetic energy of a particle impinging on 
the packing is fully dissipated by the fluid. Thus, for cohesionless particles, the dynamics of the 
deposition process, rather than the stability of the static packing, defines the random loose packing 
limit. We have made direct measurements of the interparticle friction in the fluid, and present an 
experimental measurement of the loose packing volume fraction, (J)rlp, as a function of the friction 
coefficient /is- 



INTRODUCTION 



The most elementary characteristic of a disordered 
sphere packing is the fraction, 0, of the total volume 
occupied by particles. Stable packings of cohesionless, 
frictional, spheres exist over a broad range of volume 
fractions [jj-y. The term "random close packed" (RCP) 
refers to the upper bound (J)rcp on the volume frac- 
tion at which a packing of identical spheres can be pre- 
pared without introducing crystalline ordering. Packings 
of (j) ~ 0.64, are consistently achieved in experiments 
and simulations; this number is insensitive to variations 
in interparticle forces and to the compaction protocol, 
however, questions remain as to whether there exists a 
tight upper bound [5| . The robustness of (J)rcp has mo- 
tivated attempts to understand RCP in purely geomet- 
ric terms [6| or in terms of the statistical mechanics of 
hard spheres or soft spheres at zero temper ature[7H9|. In 
this article we experimentally explore the lower bound 
on volume fractions of mechanically stable packings of 
frictional, noncohesive, identical, hard spheres. Specifi- 
cally, our questions are: does a loose packing limit exist 
for frictional but cohesionless spheres? Will this loose 
packing limit depend on the properties of the particle 
and the preparation protocol or, like the RCP limit, will 
it be robust to changes in these variables, and possibly 
admit descriptions in terms of the statistical mechanics 
of hard spheres |1Q|? 

In the first systematic study of loose packings, On- 
oda and Liniger [11] sedimented glass spheres in fluids 
of varying densities p/, approaching the density ps of 
the sphere. They found that the packing fraction ap- 
proached an asymptotic "random loose packed" (RLP) 
value, (l)RLP = 0.555, in the limit of vanishing gravita- 
tional acceleration in the fluid, gf = g {1 — pf/ps) -^ 0. 
However, the limit gf ^ conflates two different physical 
effects, both of which may plausibly lead to lower volume 
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fractions. The first effect involves the dynamics of assem- 
bling the packing: as ^/ -^ 0, falling spheres reach the 
packing with less inertia to explore the surface and re- 
arrange their neighbors. This can trap the particles in 
higher-energy, fluffier packings. A second, distinct, effect 
concerns the statics of the structure: as neutral buoy- 
ancy is approached, more fragile packings may become 
stable since the gravitational load borne by the packing 
vanishes relative to weak cohesive forces. Indeed, sim- 
ulations by Dong et al. [l2| argue that attractive van 
der W aals forces are impo rtant in stabilizing the pack- 
ings of lOnoda and Linigen at small gf . Arbitrarily low 
packing fractions can be attained when attractive inter- 
particle forces are dominant [13| , which calls into question 
the ability to experimentally access an RLP limit for co- 
hesionless spheres. 

A key goal of our experiments is to peel apart the dis- 
tinct effects of fluid properties on the statics and dynam- 
ics. Independently controlling the viscosity and density 
of the fluid allows us to test whether a unique RLP limit 
is reached as gentle deposition conditions are approached 
along arbitrary directions in the density-viscosity plane. 
Approaching the limit of gentle deposition by increasing 
viscosity while keeping the gravitational load finite al- 
lows us to avoid the cohesive regime and test whether a 
well-defined (J)rlp exists for noncohesive spheres. 

In the absence of friction, the RLP and the RCP lim- 
its are believed to coincide[14]. However, as discussed 
above, stable packings with (j) < (J)rcp are common, with 
the packing fraction showing some material-dependence 
[l|, Q. The relevant material property has been conjec- 
tured to be surface roughness 0, [ill and is experimen- 
tally found to correlate with angle of repose [15|. Thus, 
the cause of this variability is generally modelled as a fric- 
tion coefficient. Simulations with friction find the RLP 
limit to be a systematically decreasing function of friction 
[l6-l9] albeit with unexpectedly large values {p ~ 1.) 
needed to reproduce values seen experimentally [ll[ Il5j . 
In this work we directly explore (J)rlp as a function of 
measured friction coefficient with noncohesive spheres. 
We emphasize that particle contacts can exert normal 
and tangential forces, but cannot support tension. 

Another goal of our experiments, complementary to 



previous work, is to produce packings by the sequen- 
tial addition of particles, rath er than by coll ective proce- 
dures. Recent experiments bv I Jerkins et all have studied 
sett ling following brief pulses of flow in a liquid-fluidized 
bed[l5| arriving at volume fractions similar to those pro- 
duced by collective sedimentation [11|. Simulations of 
sedimented packings ^6l[l3 have studied volume fraction 
as a function of particle properties such as friction and in- 
elasticity. Other simulations 1 14 Il8| generate disordered 
packings by collectively relaxing particle conflgurations 
as volume, pressure, temperature or particle interactions 
are smoothly varied. However, packings created by se- 
quential deposition, in which a particle comes to rest at 
the flrst mechanically stable location that it encounters, 
may lead to different packings than those obtained by 
collective preparation protocols. 

In this article we present data on loose sphere packings 
prepared by the sequential sedimentation of frictional, 
non-cohesive spheres large enough to eliminate the influ- 
ence of van der Waals and other attractive forces. By us- 
ing fluids of varying density and viscosity, we identify the 
parameter in the deposition process that controls passage 
to a putative RLP limit. Rather than ^j, this parame- 
ter is the Stokes number, St = (2/9) ps Vr/r], where r is 
the radius of the sphere, V its velocity and 77 is the dy- 
namic fluid viscosity. We also vary the friction between 
spheres, both by varying the material and by increasing 
surface roughness via controlled etching. We flnd that 
the packing fraction in the loose packing limit is a func- 
tion of interparticle friction, the values of which are quite 
high, a result in qualitative agreement and simulational 
flndings. 



II. EXPERIMENTAL METHODS 

We prepared mechanically stable packings of monodis- 
perse spheres immersed in fluids, in an hour-glass shaped 
apparatus (Fig. [T]) placed on a vibration-isolation table. 
Using a variety of fluids (Table H]) and spheres (Table [III) 
we formed packings under a wide range of viscosity and 
buoyancy conditions allowing distinctions to be drawn 
between the relative merits of different parameters con- 
trolling approach to the RLP limit. 

Packings were prepared by inverting the hour-glass 
shaped cell and allowing particles to settle through the 
fluid under gravity. The hour-glass geometry consists of 
two conical sections with a cone angle of 60 ° and base 
diameter of 24 sphere-diameters {d = 2r) connected by 
a cylindrical neck of 4.2 d in diameter, which is only as 
wide as necessary to avoid jamming by arch formation in 
the neck. The narrow neck allows the passage of only a 
few particles at a time. (We have also deposited spheres 
singly by a mechanical dropper, with very similar results 
to those reported here.) The packing grows as a conical 
pile at the angle of repose (~ 23 — 26°), which is much 
smaller than the cone angle of the container, thus elim- 
inating any empty pockets near the walls. The conical 




FIG. 1: Hour-glass shaped sedimentation apparatus. Top 
right: Image of the topography of the top surface of the par- 
ticles. Since the volume is estimated from a single projection, 
there is a small systematic positive bias of ^0 < +0.002 in 
measuring the volume of the packing. 
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FIG. 2: Histograms of roughness C, of spherical particles as 
measured with a DektakS profilometer. Inset (a): Select pro- 
filometer traces C as a function of distance along the surface x 
(the gross curvature of the sphere has been subtracted) . Im- 
ages of spherical particles: (A) steel, (B) smooth acrylic, (C) 
PTFE, (D) solvent-etched acrylic, (E) aluminum. 



cell was chosen to minimize the weight supported by the 
sloping walls of the container. No crystalline order was 
observed near the bottom or side walls. 

Data for the volume fraction of the packing are taken 
when the top surface of the packing just enters the neck 
of the cell. The total volume of the packing is the volume 
of the cone plus the small contribution from the spheres 



TABLE I: Fluids used in hour-glass sedimentation experiments. The symbols match those used in all figures. Mixture ratios 
reflect contents before degassing. Density and viscosity were measured after degassing. *lit.|20l]. ^our measurement, "'"lit. at 
20C[2i. §MSDS. 



fluid 


density /{gem ^) 


viscosity/ (mPa s) 


symbol 


n-pentane (Fisher Sci) 


0.626* 


0.23* 


M 


n-heptane (Fisher Sci) 


0.684* 


0.41* 


> 


water (Millipore) 


1.00* 


1.03* 


V 


n-dodecane (Fisher Sci) 


0.75 


1.53* 


A 


^-^71/29 n-dodecane/light mineral oil mixture 


0.78* 


3.02* 


o 


^^50/50 n-dodecane/light mineral oil mixture 


0.80* 


5.96* 




^-^48/52 n-dodecane/light mineral oil mixture 


0.82* 


11.5* 




light mineral oil (Fisher Sci) 


0.83§ 


46.0* 





high temperature silicone oil (Acros Organics) 


1.05§ 


117.* 




O.OIM NaCl in --70/30 propylene glycol/glycerol 


1.09* 


125.* 




heavy mineral oil (Fisher Sci) 


0.83§ 


157.* 




Fluka 08577 Density and Viscosity Standard 


0.87§ 


1270.§ 


T 



TABLE IL Properties of sets of spheres. We measured particle diameter using a technique similar to Scott'sjj], measuring the 
length of ^ 100 spheres in a groove. Density was calculated from this diameter and the weight of these spheres. To quantify 
polydispersity and sphericity, we measured the diameter of individual spheres with a machinist's micrometer accurate to 2.5/im, 
along flve or more directions. "Asphericity" is the relative deviation from sphericity calculated as the standard deviation of 
these diameter measurements relative to the mean. "Polydispersity" is the standard deviation of the average diameters of a set 
of 20 spheres relative to the mean. "RMS roughness" gives the root-mean-square deviation from sphericity of profilometer traces 
taken of the sphere's surface(see Fig. [2]). fis and a^^ are the mean and width of the distribution of static friction coefficient of 
gently contacting spheres in fluid (see Fig. [5|). 



sphere material 


diameter/cm 


density /{gem ^) 


polydispersity 


asphericity 


RMS roughness//im 


fIs 


C^/x, 


PTFE 


0.3205(4) 


2.1389(9) 


0.21 ±0.06% 


0.14 ±0.06% 


1.1 ±0.6 


0.540 ± 0.003 


0.10 


aluminum 


0.3191(3) 


2.775(3) 


< 0.04% 


< 0.04% 


0.32 ±0.14 


0.62 ± 0.06 


0.16 


steel 


0.3179(4) 


7.774(7) 


0.14 ±0.05% 


< 0.04% 


0.10 ±0.02 


0.66 ±0.14 


0.15 


smooth acrylic 


0.3174(3) 


1.1800(9) 


0.15 ±0.06% 


0.06 ± 0.04% 


0.7 ±0.3 


0.88 ± 0.03 


0.10 


etched acrylic 


0.3092(3) 


1.1741(9) 


0.16 ±0.07% 


0.08 ± 0.05% 


2.6 ±0.1 


0.96 ± 0.03 


0.10 



that are in the neck. The latter is determined from an 
image of the topography of the top surface of the particles 
(inset Fig.[T]). All volume fractions reported in this article 
are subject to the same systematic error in the range 
S(j) ~ 0.000-0.002 due to uncertainties in these volumes 
and in the volume of the hour glass. 

The particles used in our experiments are com- 
mercially available spheres of acrylic (PMMA), teflon 
(PTFE), steel, and aluminum with nominal diameter 
d ^ 3.18mm(see Table HH Fig. [2j). We also use acrylic 
spheres which were roughened by timed etches in acetone. 
All sets of spheres are very monodisperse and highly 
spherical with, at worst a deviation from sphericity of 
~ 10~^(i, a surface roughness of the same order, and a 
polydispersity of double this magnitude. These spheres 
are much larger than those used by Onoda and Liniger 
(0.25 ± 0.02 mm, glass). The advantages of using large 
spheres are evident: attractive van der Waals forces are 
negligibly small compared to other forces in the problem 
and the particles are well characterized and of extremely 
high sphericity and monodispersity. Thus these experi- 



ments represent a better approximation to the idealized 
packing of hard, monodisperse spheres than previous ex- 
periments. 

Apart from van der Waals attr action [l2|, experimen- 
tal results on the low packing fraction limit are sensi- 
tive to other attractive forces of electrostatic and cap- 
illary origin. In setting up the apparatus, great care 
was taken to degas fluids and to introduce the fluid to 
the spheres slowly to avoid entraining air bubbles which 
form attractive bridges between between poorly wetted 
surfaces, especially rough ones. Pumping a vacuum on 
heated, stirred fluids for hours was often insufficient to 
avoid the appearance of attractive forces between spheres 
in non- wetting fluids, a phenomenon that has recently 
been associated with the existence of long-lived [2l|, |22| 
nanobubbles capable of exerting forces comparable to 
gravity for PMMA spheres in our more closely density 
matched fluids [23|. To avoid this phenomenon, we have 
used well- wetting fluids when possible and avoided close 
density matches in poorly- wetting fluids so that the con- 
tribution of attractive forces is negligible. Where charg- 



0.59 



0.58 



-^ 0.57 



eO.56 



0.55 - 



0.54 



0.53 




'1^ 



0.59 



(c) 






J ^ ^ ^ ^ ^ ^ 0.53 



10"^ . 10"' 



10" 



R-^ =^rfl{2r^ PfPsQf) 



(b) 



« 'V 



10-^ 



0.59 
0.57 

4 

0.55 
0.53 



10- 






lO'^ 



12 4 8 

Reduced Gravity, (g^/g)"^ 



FIG. 3: Packing fraction (j) vs. the dimensionless parameters S^ R, and Qf/g. The data points represent individual packings; 
the spread in the data is much larger than the random error on each data point. For a particle falling at low Reynolds number, 
S is the Stokes number and R is the Reynolds number itself, gf is the buoyancy-reduced gravitational acceleration felt by a 
particle in the fluid. (Onoda and Liniger's A^ = 5'//pJJj].) In all three graphs the solid line connects data for smooth acrylic 
spheres and the dashed line is for steel spheres. 



ing effects were suspected, we repeated measurements 
with salts added to screen coulomb interactions. 



III. RESULTS AND DISCUSSION 

The problem of a sphere falling in the fluid involves 
five dimensionful parameters: r = d/2 and ps, the radius 
and density of the spheres; pf and r^, the density and dy- 
namic viscosity of the fluid; and gf = g {1 — pf/ps)^ the 
buoyancy-reduced gravitational acceleration in the fluid. 
Apart from gf/g (as suggested by Ref. [ll|), other perti- 
nent dimensionless groups are the Reynolds number Re = 
2pfVr/r] and the Stokes number St = {2/9) psV r/r]^ 
where V is the velocity of the particles as they approach 
the packing. In reducing our data we use for V^ the ter- 
minal velocity of the particle under Stokes drag, leading 
to the dimensionless parameters R = (2/9) r^ ps pf gf/rf 
and 5* = {4/81) r^ p^ gf/rj'^. In the limit of gentle depo- 
sition, R = Re^ S = St. S can be interpreted either as a 
dimensionless damping length, or as the ratio of kinetic 
energy to the potential energy, mass x gf d, which quan- 
tifies the degree to which a falling sphere can rearrange 
the packing. 

We display in Fig. [3] the major result of this paper. 
For spheres of a given material, the packing fraction 
decreases as sedimentation is done more gently. When 
plotted against the dimensionless group 5, the volume 
fraction approaches an asymptotic limiting low volume 
fraction which we may interpret as (J)rlP' We note the 



the limiting (J)rlp is different for spheres of different ma- 
terials. The asymptotic limit is directly available from 
the data, and unlike in previous measurements need not 
be obtained by extrapolation (for which there is no re- 
liable theoretical guideline). Each data point is taken 
with a different liquid, and not with a chemical series or 
a dilution series; the smooth approach to this limit thus 
suggests that the macroscopic parameters of the fluid and 
sphere are sufficient to fully characterize the preparation 
process, and that any microscopic interparticle interac- 
tions mediated by the fluid have been successfully sup- 
pressed. 

Furthermore, Fig. [31(b) shows that the degree of den- 
sity matching, quantified hy gf/g = 1 — pf/ps^ is not 
an appropriate parameter to define the RLP limit of this 
packing protocol. Low packing fractions can be achieved 
without density matching, by sedimenting in sufficiently 
viscous fluids. For a fluid of a given viscosity, RLP can 
of course be approached by varying gf/g^ as the work of 
Onoda and Liniger suggests. The role of gf/g is however, 
not in reducing the static load on the packing structure, 
but purely in slowing down the dynamics of the packing 
process. 

Fig. [31(c) shows that the volume fraction (j) also 
smoothly approaches (I)rlp when plotted against the di- 
mensionless group R. However, unlike the plot of (/) vs. 
5, the asymptotic limit is attained at different values of 
R for different materials, leading us to prefer S as the 
best candidate for the relevant control parameter. 

Additionally, previous measurements of collisions of 
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FIG. 4: The main plot shows the packing fraction, 0, of 
acryhc spheres (left vertical axis, symbols match Fig. [3]) ver- 
su s Stokes numbe r, St. We compare the trend with the data 
of lGondret et a/.l (from Fig. 6, Ref. [23]), plotted here in dark 
circles against the right vertical axis) for the 5't-dependence of 
the restitution coefficient e, scaled by its maximal value Cmax- 
The data are for collisions of a teflon sphere in a fluid, how- 
ever, data for different materials collapse on the same curve. 
To make the comparison with Ref. 24] , we inferred the value 
of St in our experiments from measured values f Table HI)) and 
a standard drag curve 25]. This correction is shown in the 
inset; as expected the correction is small when St is small. 



FIG. 5: Limiting low packing fraction (I)rlp plotted versus /is, 
the mean static friction coefficient, /is is measured by a new 
technique where the friction between spheres immersed in a 
fluid is ascertained from the maximum angle at which pairs 
of spheres in a vertical, circular track maintain contact due 
to static friction. The dashed error bars indicate the width in 
the distribution of /j^s , and the solid error bars are uncertainty 
in the mean of /Js- The shaded area represents the region of 
stable, disordered packings. The lower bound (dashed curve) 
is a guide to the eye. The inset shows a drawing of the ge- 
ometry used in the friction measurement (scale: bead radius 
r = 1.59 mm). 



fluid- immersed spheres on surfaces [24j also show that the 
Stokes number St - and not Reynolds number Re - is the 
relevant dimensionless parameter that defines the onset 
of bounce-free collisions. In addition to supporting the 
physical picture that the RLP limit corresponds to condi- 
tions where sequentially added spheres settle in the first 
mechanically stable location that they encounter, there 
is good quantitative agreement between the scale of St 
where bouncing ceases, and where (J)rlp is attained. This 
correspondence is shown in Fig. HJ The Stokes number 
has also been identified as the parameter controlling the 
behavior of avalanches in fluid- immersed piles [26|]. 

We now return to the observation that the two curves 
in Fig. [3] corresponding to steel and acrylic spheres ap- 
proach different values oi(j)RLP of 0.551 and 0.540 respec- 
tively. With the mass density of the materials already 
accounted for, the only relevant differences are those of 
the contact mechanics of the spheres. In particular, the 
effective coefficient of static friction /i^ between spheres 
is different for these materials. Thus, unlike the RCP 
limit, the RLP limit is not a purely geometric problem, 
but involves the mechanics of interactions. 

In order to more fully explore this observation, we 
have prepared packings of five different materials (Ta- 
ble ini) in the St ^ {) limit. The coefficient of static 
friction, /ig is affected both by material as well as by 
surface topography; indeed we find that (J)rlp does not 
show a simple trend with surface roughness. To directly 
probe the material property of interest, we have devised 
a method to measure /ig for sphere-on-sphere contacts 



between spheres immersed in a liquid. The schematic di- 
agram inset in Fig. [5] shows the geometry of the setup: 
two spheres sit under gravity in the shallow v-groove of 
a track in the vertical plane. As the track rotates slowly, 
static friction between beads prevents them from rolling 
so they move with the track until they reach a maximal 
angle Omax- At Omax tfie tangential force between the 
spheres exceeds the maximum value allowed by the fi- 
nite coefficient of static friction jig and the beads roll to 
a lower an^le Omin < Omax- l-^s can be then calculated 
from 6>macc|23- I^ Fig. [5] we show (J)rlp decreases mono- 
tonically as the measured friction coefficient increases. 
The dependence on friction coefficient, 4)rlp{ij^s)-, is con- 
sistent with (j)RLP -^ (j^RCP as the coefficient of static 
friction, /ig ^ 0. 

Our measurements of /ig are made under conditions 
similar to those in our packing experiments, under the 
same small loads and fluid environments, and with 
sphere-on-sphere contacts that allow for both sliding and 
rolling. We are not aware of any other measurements of 
interparticle friction in this regime. The values of /ig we 
observe are larger than those from our everyday expe- 
rience at larger normal loads. Our data thus suggest a 
resolution of the puzzle that friction coefficients in sim- 
ulated contact mechanics models |16Ml8| were thought to 
be surprisingly large in order to achieve volume fractions 
as low as seen in experiments. 

Finally, we turn from the RLP limit to the large Stokes 
regime. Surprisingly, we see a plateau in the value of (/), 
well separated from the low Stokes limit. The plateau 
(j) does not scale simply with friction coefficient /i^, and 



presumably also involves other particle properties such 
as inelasticity. We speculate that this plateau value may 
be related to the "critical state" of soil mechanics. The 
transition from this plateau to (J)rcp is clearly of interest, 
but is not easily explored by merely tuning S with fluid 
parameters. 



IV. CONCLUSION 

The strengths of our experiments are (i) the extremely 
well-controlled sphericity and monodispersity of the par- 
ticles, (ii) detailed characterizations of relevant particles 
properties such as surface roughness and of friction co- 
efficient under deposition conditions, (iii), broad cover- 
age of fluid parameters, and (iv) employing large enough 
particles to be well outside the influence of any attrac- 
tive interactions. Thus our experiments provide the best 
available approximation to the idealized problem of the 
packing of monodisperse frictional spheres. Despite our 
choice of experimental geometry, the relatively small sys- 
tem size in our experiments may introduce wall-effects. It 
is clear, however, that the values of (I)rlp obtained by the 
sequential deposition in our experiment are comparable 
with those obtained by collective packing schemes [ill [l5j , 
and therefore our findings should also be applicable to 
those packing protocols. 

Previous experiments have shown that looser packings 



may be prepared if sedimentation is done more gently. 
Our results strengthen this intuitive expectation in three 
significant directions. The first is that we arrive at a 
sharp definition of "gentle" deposition: this limit is gov- 
erned by the Stokes number, St. The second conclusion 
is that the RLP limit is achieved at a nonzero threshold 
value of St., below which particles entering the packing do 
not have the ability to explore the landscape of deposited 
particles, or to rearrange it. Finally, since we eliminate 
the effect of attractive forces by packing at finite values of 
^/, we establish the existence of a (I)rlp for cohesionless 
spheres. 

We also provide the first direct experimental study of 
the friction-dependence of (I)rlp by measuring the fric- 
tion coefficients between particles at small normal loads. 
Friction stabilizes packings at volume fractions consider- 
ably below (j)RCP' The evolution of a packing from the 
RLP boundary to the RCP boundary, and the structure 
and mechanical properties of the intermediate states re- 
main largely unexplored. 
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